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We report the first study of the glueball properties at finite temperatures below T c using SU(3) 
anisotropic lattice QCD with (3 = 6.25, the renormalized anisotropy £ = a s /a t = 4 and 20 3 x Nt 
(-/Vt=35,36,37,38,40,43,45,50,72) at the quenched level. From the temporal correlation analysis with 
the smearing method, about 20 % mass reduction is observed for the lowest scalar glueball as 
m G (T) = 1250 ± 50MeV for 0.8T C < T < T c in comparison with m G ~ 1500~1700MeV at T ~ 0. 
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Finite temperature QCD, including the quark gluon 
plasma (QGP) physics, is one of the most interesting 
subjects in the quark hadron physics jy, [|, g]. At high 
temperature, in accordance with the asymptotic freedom 
of QCD, the strong interaction among quarks and gluons 
is expected to be reduced, and there would occur the 
deconfmement and/or chiral phase transition 

For the study of finite temperature QCD, the lattice 
QCD Monte Carlo simulation provides a reliable method 
directly based on QCD. For instance, SU(3) lattice QCD 
simulations at the quenched level show a weak first-order 
deconfmement phase transition at the critical tempera- 
ture T c ~ 260MeV {§, and full SU(3) QCD simulations 
show a chiral phase transition at T c = 173(8)MeV for 
N f = 2 and 154(8)MeV for N f = 3 in the chiral limit [§. 
Above T C) most of the nonperturbative properties such 
as color confinement and spontaneous chiral-symmetry 
breaking disappear, and quarks and gluons are liberated. 

Even below T c , there are many model predictions on 
the change of the hadron properties^], ^, the mass and 
the size, due to the change in the inter-quark potential 
H [| and the partial chiral restoration. As a pre-critical 
phenomenon of the QCD phase transition, the possible 
hadron mass shift at the finite temperature or in the fi- 
nite density is now one of the most interesting subjects 
in hadron and QGP physics. For instance, the CERES 
data with the ultra-relativistic heavy-ion collision experi- 
ment may indicate the p-meson mass shift JToj ] , and many 
theoretical studies [0 have been done to explain this ex- 
periment. 

Nevertheless, lattice QCD studies for thermal proper- 
ties of hadrons are still inadequate at present because of 
the difficulty in measuring the hadronic two-point corre- 
lators on the lattice at finite temperature. For instance, 
on the screening-mass measurement 113] , this difficulty 



is due to the mixture of the large Matsubara frequencies 
in addition to the absence of technical prescriptions as 
the smearing method. On the other hand, on the pole- 
mass measurement, while it is free from the mixture of 
the Matsubara frequencies, another difficulty arises from 
the shrink of the physical temporal size 1/T at high tem- 
perature. In fact, the pole-mass measurements have to 
be performed within the limited distance shorter than 
1/(2T), and such a limitation corresponds to Nt — 4 ~ 8 
near T c in the ordinary isotropic lattice QCD ||. 

To avoid this severe limitation on the temporal size, 
we adopt an anisotropic lattice where the temporal lattice 
spacing at is smaller than the spatial one a s ]8|, [l~4|, [l5| . 
We can thus efficiently use a large number of the temporal 
lattice points as Nt ~ 32 even near T c , while the physical 
temporal size is kept fixed 1/T = N t at- In this way, the 
number of available temporal data is largely increased, 
and accurate pole-mass measurements from the temporal 
correlation become possible |l5| . 

In this paper, we study the glueball at finite temper- 
ature from the temporal correlation analysis. We use 
SU(3) anisotropic lattice QCD at the quenched level, as 
a necessary first step before attempting to include the 
effects of dynamical quarks in the future. Even with- 
out dynamical quarks, quenched QCD can reproduce well 
various masses of hadrons, mesons and baryons, and im- 
portant nonperturbative quantities such as the confining 
force and the chiral condensate. In quenched QCD, un- 
like full QCD, the elementary excitations are only glue- 
balls in the confinement phase below T c ~ 260MeV. 
At zero temperature, the lightest physical excitation 
is a scalar glueball with J PC = ++ with the mass 
m G ~ 1500~1700MeV |§ [H| [HJ |J|, which is expected 
to dominate the thermodynamical properties below T c . 

We consider the glueball correlator (ll| [l6], [l7], [l^, 



2 



g |g in SU(3) lattice QCD as G(t) = (O(t)O(O)}, 
0(f) = 0{t) - (O), 0(t) = J20(t,x). The summa- 

X 

tion over x physically means the zero-momentum pro- 
jection. The glueball operator 0(t, x) is to be properly 
taken so as to reproduce its quantum number J PC in 
the continuum limit. For instance, the simplest com- 
position for the scalar glueball is given as 0(t,x) = 
ReTr {Pi 2 (t, x) + P 23 (t, x) + P 3 i(t, x)}, where P^(t, x) e 
SU(3) denotes the plaquette operator. With the spec- 
tral representation, G(t) is expressed as G(t)/G(0) — 
Ecne- E "\ c = |(n|O|0)| 2 /G(0), 0(0) = £|(n|O|0)| 2 , 
where E n denotes the energy of the n-th excited state 
\n). Here, |0) denotes the vacuum, and |1) denotes the 
ground-state glueball. Note that c n is a non-negative 
number with ^c„ = 1. On a fine lattice with the spac- 
ing a, the simple plaquette operator Pij(t, x) has a small 
overlap with the glueball ground state \G) = |1), and 
the extracted mass looks heavier owing to the excited- 
state contamination. This small overlap problem origi- 
nates from the fact that 0(t, x) has a smaller "size" of 
0(a) than the physical peculiar size of the glueball. This 
problem becomes severer as a — > 0. We thus have to im- 
prove 0(t,x) so as to have almost the same size as the 
physical size of the glueball. 

One of the systematic ways to achieve this is the smear- 
ing method |2(], ^2j. The smearing method is ex- 
pressed as the iterative replacement of the original spatial 
link variables Ui(s) by the associated fat link variables, 
Ui(s) £ SU(3) C , which is defined so as to maximize 



ReTr 



Ul(s)[aUi(s)+ ]T U ±j (s)Ui(s ±j)U{ j (s + i) 



(1) 



where U-^(s) = f/I(s — /t), and a is a real parameter. 
Here, the summation is taken only over the spatial di- 
rection to avoid the nonlocal temporal extension. Note 
that Ui(s) holds the same gauge transformation proper- 
ties with Ui(s). We refer to the fat link defined in Eq. (|l|) 



as the first fat link U^\s). The n-th fat link E^" ; (s) is 



(«), 



defined iteratively as U- 
r(l) = 77 



(n) 



(s) 



1 (s) staring from 



U^ 1 ' = Ui(s) [j22j. For the physically extended glueball 

operator, we use the n-th smeared operator, the plaque- 

(„) 

tte operator constructed with U i (s). 

The smeared operator physically corresponds to an ex- 
tended composite operator with the original field variable 
as Ufj,(s). We consider the size of the n-th smeared op- 
erator in terms of the original field variable. Using the 
linearization on the gluon field, we obtain the diffusion 
equation as (fLol p3| 



d 

t-K(x, n) = DAK(x, n), D 
on 



s 



(2) 



for the distribution K(x, n) of the gluon field in the n- 
th smeared plaquette, in the case of the small spatial 



TABLE I: The lattice QCD result for the lowest scalar glueball 
mass at finite temperature. The temporal lattice size Nt, 
the corresponding temperature T, the lowest scalar glueball 
mass mc(T), the maximal ground-state overlap C ,max ; fully 
correlated \ 2 /Ndf, the smearing number N BU1T , the number 
of gauge configurations iV con f and the rough estimate of the 
glueball size p are listed. The most suitable smearing number 
-/V smr is determined with the maximum ground-state overlap 
condition. 



Nt 


T[MeV] 


m G [MeV] 




XV^DF 


smr 


N con f 


pM 


72 


130 


1450(40) 


0.93(2) 


1.43 


39 


5541 


0.42 


50 


187 


1410(46) 


0.92(3) 


0.34 


41 


5168 


0.44 


45 


208 


1456(34) 


0.96(1) 


0.72 


40 


5929 


0.43 


43 


218 


1323(39) 


0.89(2) 


0.90 


43 


8693 


0.45 


40 


234 


1260(45) 


0.84(3) 


0.75 


42 


7420 


0.44 


38 


246 


1221(35) 


0.85(2) 


0.12 


40 


8736 


0.43 


37 


253 


1273(32) 


0.88(2) 


1.61 


38 


8633 


0.42 


36 


260 


1208(35) 


0.84(2) 


1.34 


39 


8603 


0.42 


35 


268 


1188(34) 


0.84(2) 


1.80 


40 


8462 


0.43 



lattice spacing a s . The n-th smeared plaquette located 
at the origin x = physically correspond to the Gaussian 
extended operator with the distribution as E3[ 



K(x, n) = 



1 



(7^2)3/: 



■ exp 



^21 



(3) 



where p represents the characteristic size of the Gaussian 
distribution, and is defined as 



2a. 



(4) 



Thus, the smearing method, which is introduced to carry 
out the accurate mass measurement by maximizing the 
ground-state overlap, can be also used to give a rough 
estimate of the physical glueball size. In fact, once we 
obtain the maximum overlap with some n and a, the 
glueball size is roughly estimated with Eq. (|^). 

We use the SU(3) anisotropic lattice plaquette action 



S e = W— E ReTr(l-^(s)) 



(•5) 



7G ReTr(l-P l4 (s)) 



s , i < 3 



with the plaquette operator P„ w (s) £ SU(3) in the (//, v)- 
plane. The lattice parameter is fixed as (3 = 2N c /g 2 = 
6.25, and the bare anisotropy parameter is taken as 7g = 
3.2552 so as to reproduce the renormalized anisotropy 
£ = a s /a t = 4 ]l3| . These parameters produce the spa- 
tial lattice spacing as aj 1 = 2.341(16) GeV (a s ~ 0.084 
fm), and the temporal one as a^ 1 — 9.365(66) GeV 
(at ~ 0.021 fm). Here, the scale unit is determined 
by adjusting the string tension as y/cr = 440MeV from 
the on-axis data of the static inter-quark potential. The 
pseudo-heat-bath algorithm is used to update the gauge 



3 



1 



field configurations on the lattice of the sizes 20 3 x N t , 
with N t = 35, 36, 37, 38, 40, 43, 45, 50, 72 as listed in Table 
For each temperature, we pick up gauge field config- 
urations every 100 sweeps for measurements, after skip- 
ping more than 20,000 sweeps of the thermalization. The 
numbers of gauge configurations used in our calculations 
are summarized in Table [j]. 

For completeness, we give an estimate of the crit- 
ical temperature T c . To this end, we analyze 
the scattering plot of the Polyakov loop P{x) = 
TilU^x, 0) • • • U±{x, N t — 1)} at each gauge field config- 
uration. From this analysis, the Z 3 symmetry holds at 
N t — 35, and the system is found to be in the confinement 
phase. On the other hand, the Z3 symmetry is broken 
at Nt — 34, which indicates the deconfincmcnt phase. 
Hence, we estimate T c ~ 270MeV, which is consistent 
with the previous studies j|, |) . 

We present the numerical results in SU(3) anisotropic 
lattice QCD at the quenched level. To enhance 
the ground-state contribution, we adopt the smearing 
method with the smearing parameter a — 2.1, which we 
find one of the most suitable values from the numerical 
tests with various a. The statistical errors are estimated 
with the jackknife analysis ||l9|| . 

In Fig. 0(a), we show a scalar glueball correlator 
G(t)/G(0) at a low temperature T = 130MeV for the 
smearing number 2V SBlr = 40, where most of the lattice 
QCD data are well fitted by a single hyperbolic cosine, 
denoted by the solid curve, as 



(a) 



G(t)/G(0) = C(e' mGtat + e -™G(N t ~t)a t y 



(6) 



This indicates the achievement of the ground-state en- 
hancement owing to the smearing method, and then the 
excited-state contamination is almost removed. 

In general, G(t)/G(0) is expressed as a weighted sum 
of hyperbolic cosines with non-negative weights, and 
G(t)/G(0) decreases more rapidly than Eq. (||) near t = 
due to excited-state contributions. Hence, C should sat- 
isfy the inequality C < (1 + e -™-oa t N t yi _ x In thc 

ground-state dominant case, G(t)/G(0) can be well ap- 
proximated by a single hyperbolic cosine, and C ~ 1 is 
realized. We refer to C as the ground-state overlap. 

From Fig. 0(a), we find C ~ 1 and m G ~ 1450MeV 
for the lowest scalar glueball mass at a low temperature. 
This seems consistent with m G ~ 1500~1700MeV at T ~ 
|1|. 

In Fig. 0(a), we show a scalar glueball correlator 
G{t)/G(0) at a high temperature T = 253MeV for the 
smearing number iV smr = 40. Owing to a suitable smear- 
ing, most of the lattice QCD data are well fitted by a 
single hyperbolic cosine denoted by the solid curve. 

Each best fit analysis is performed in the interval 
[tmin> t max ], which is determined from the flat region 
[£min,^max — 1] appeared in the corresponding "effective 
mass" plot shown in Figs. [11(b) and 0(b). The effective 
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FIG. 1: (a) The scalar glueball correlator G(t)/G(0) for 
jV smr = 40 at a low temperature T = 130MeV. (b) The cor- 
responding effective mass plot. The statistical errors are es- 
timated with the jackknife analysis. The solid line denotes 
the best single hyperbolic cosine fit to the lattice data in the 
interval [t mm , tmax] indicated by the two vertical dashed lines. 
The dashed and dotted curves are the best hyperbolic cosine 
curves for the modified fit range with t m j n + 1 and + 2, 
respectively. The closeness of the three curves means small 
fit-range dependence. 



mass m e g(t) is a solution of 

G(t + 1) cosh (m eS {t)a t (t + 1 - N t /2)) 



G(t) 



cosh (m e s(t)at(t — Nt/2)) . 



(7) 



for agiven G(t + 1)/G(t) at each fixed t @. In Figs. 
and we show also the results of further two fits in the 
modified interval as [imin + Mmad and [t min + 2, t max ] by 
dashed line and dotted line, respectively. The closeness 
of the three curves suggests small fit-range dependence. 

In the most suitable smearing iV smr , the ground-state 
overlap C is maximized and the mass uiq is minimized, 
which indicates the achievement of the ground-state en- 
hancement. (For extremely large N sxal , the operator 
size exceeds the physical glueball size, resulting in the 
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FIG. 2: (a) The scalar giueball correlator G(t)/G(0) for 
jV smr = 40 at a high temperature T = 253MeV. (b) The 
corresponding effective mass plot. The statistical errors are 
estimated with the jackknife analysis. The solid line denotes 
the best single hyperbolic cosine fit to the lattice data in the 
interval [t m in, imax] indicated by the two vertical dashed lines. 
The dashed and dotted curves are the best hyperbolic cosine 
curves for the modified fit range with t m j n + 1 and t m i n + 2, 
respectively. The closeness of the three curves means small 
fit-range dependence. 



decrease of the overlap C.) In practical calculations, 
the maximum overlap and the mass minimization are 
achieved at almost the same N sml , and both of these two 
conditions would work as an indication of the maximal 
ground-state enhancement. Here, we take the maximum 
ground-state overlap condition as C ~ 1. (The mass 
minimization condition leads to almost the same glue- 
ball mass [p3|.) 

From the analysis at various temperatures, we plot the 
lattice QCD result for the lowest scalar giueball mass 
ma(T) against temperature T in Fig. ||. We observe, 
in Fig. ||, about 20 % mass reduction or a few hundred 
MeV mass reduction of the lowest scalar giueball near 
T c as m G (T) = 1250 ± 50MeV for 0.8T C < T < T c in 



I 



> 



I 



I 



I 



1000 - 



150 



I I L_ 

200 
7TMeV] 



250 



FIG. 3: The lowest scalar giueball mass plotted against the 
temperature T. It is obtained with the best hyperbolic cosine 
fit in the interval [i m i n ,imax] determined from the flat region 
in the effective mass plot. The vertical dotted line indicates 
T c ~ 260MeV. 



comparison with m G {T ~ 0) ~ 1500-1700MeV g§, |T|, 

g. 

We also give a rough estimate of the giueball size. To 
estimate the giueball size, we search N sml which realizes 
the maximum ground-state overlap C max . From Eq. (^) 
with this iV smr , we roughly estimate the giueball size as 
p ~ 0.4 ~ 0.45 fm both at low temperature and at high 
temperature near T c . Thus, we see that the thermal effect 
on the giueball size is rather small, which may provide 
an important information in the bag model argument of 
the QCD phase transition Jl5|, |23| . 

In Table |, we summarize the lowest scalar giueball 
mass mc{T), the ground-state overlap (7 max , fully corre- 
lated xV-^df, the corresponding smearing number iV smr , 
the number of gauge configurations N con { and the esti- 
mated giueball size p. 

Thus, the present lattice QCD calculation indicates 
that the lowest scalar giueball exhibits about 250MeV 
mass reduction near T c keeping its size. Here, we briefly 
discuss the physical consequence of this result, consider- 
ing the trigger of the QCD phase transition. In quenched 
QCD below T c , the lowest giueball is the lightest parti- 
cle, and its thermal excitation is expected to have pri- 
mary relevance at finite temperature. However, lattice 
QCD indicates toq > lGeV even near T c , and there- 
fore the thermodynamical contribution of the giueball 
seems strongly suppressed by the small statistical factor 
as e -" lG / T near T c ~ 260MeV ||. This may indi- 
cate that the thermal giueball excitation does not play 
the relevant role in the deconfinement phase transition, 
at least in quenched QCD. Then, what is the driving 
force to bring the phase transition ? In this way, our 
result brings up such an interesting new puzzle on the 
QCD phase transition. 

Several comments are in order. The first comment 
is on the closeness of our simulations to the continuum 
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limit. In Ref.[Q, the authors investigated f3 dependence 
of glueball masses at zero temperature, and estimated 
the discretization error on the scalar glueball mass to 
be less than 5 % at f3 — 6.4. According to them, the 
discretization error is estimated about 6 % at (3 = 6.25 
in the present calculation. The second comment is on 
the finite volume artifact on the scalar glueball mass. In 
Rcf. (T^| , Monte Carlo simulations were performed on the 
lattice of the physical size (1.76fm) 3 and (1.32fm) 3 at 
zero temperature to investigate the finite volume errors 
in the various glueball masses by using an improved ac- 
tion. The authors concluded that the systematic error 
in the lowest scalar glueball mass from the finite volume 
is negligible at zero temperature. Note that the finite 
volume artifact on the scalar glueball mass is essentially 
independent of the regularization method, i.e., a specific 
choice of the lattice action, as far as the discretization is 
enough fine. It follows that, the finite volume artifact of 
our results are negligible, since the physical size of our 
lattice is (1.68fm) 3 . 

To summarize, we have studied the glueball properties 
at finite temperature using SU(3) anisotropic quenched 
lattice QCD with 5,000-9,000 gauge configurations at 
each temperature. From the temporal correlation anal- 
ysis with the smearing method, we have observed about 



20% mass reduction of the lowest scalar glueball as 
m G (T) = 1250 ± 50MeV for 0.8T C < T < T c , while no 
significant ch ang e is seen for meson masses near T c in 
lattice QCD JTj. 

Finally, we comment the brief outlook. It seems in- 
teresting to investigate other glueballs such as the 2 ++ 
glueball at finite temperature to clarify whether the ther- 
mal mass reduction is peculiar to the lowest scalar glue- 
ball or universal feature in glueballs. It would be also 
interesting to analyze the spectral function of the glue- 
ball at finite temperature from its temporal correlation 
in terms of the mass and the thermal width, because 
the width broadening may provide the similar effect to 
the temporal correlator jllj] as the mass reduction. Our 
result shows that the scalar glueball mass reduction is 
about 250MeV, which is enough large, and therefore the 
thermal mass shift of the scalar glueball may become ob- 
servable in the future experiment in RHIC. 
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